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Electrons hopping on the sites of a three-dimensional pyrochlore lattice are shown to form topo- 
logically non-trivial insulating phases when the spin-orbit (SO) coupling and lattice distortions are 
present. Of 16 possible topological classes 9 are realized for various parameters in this model. Specif- 
ically, at half-filling undistorted pyrochlore lattice with SO term yields a 'pristine' strong topological 
insulator with Z2 index (1;000). At quarter filling various strong and weak topological phases are 
obtained provided that both SO coupling and uniaxial lattice distortion are present. Our analysis 
suggests that many of the non-magnetic insulating pyrochlores could be topological insulators. 
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According to recent pioneering theoretical studies [T] |2] 
all time-reversal (T) invariant (non-magnetic) band in- 
sulators in three spatial dimensions can be classified into 
16 topological classes distinguished by a four-component 
topological index {vq',viI'2V3) with = 0,1. Ordinary 
'trivial' band insulators have index (0;000) and, in gen- 
eral, possess no robust surface states. When some of 
the vs differ from zero then the insulator is said to be 
topologically non-trivial and, as a result, possesses topo- 
logically protected surface states on at least some of its 
surfaces. When — \ surface states exist on all sur- 
faces and are in addition robust with respect to weak 
non-magnetic disorder. This is referred to as a strong 
topological insulator (STI) . Strong topological insulators 
are predicted to exhibit a host of unusual phenomena 
associated with their non-trivial surface states. These 
include proximity-induced exotic superconducting state 
with Majorana fermions bound to a vortex [3], spin- 
charge separated solitonic excitations [11 15], and, in a thin 
film geometry, an unconventional excitonic state with 
fractionally charged vortices [B]. There are also interest- 
ing bulk manifestations of STI physics such the 'axion' 
electromagnetic response [71 [S] and the topologically pro- 
tected fermion modes localized along the core of a crystal 
dislocation [9J. 

Topologically nontrivial insulating phases have been 
predicted to occur |Tni HU [H] and subsequently ex- 
perimentally discovered (THl EH |TS] in several 2 and 3- 
dimensional crystalline solids. Vigorous search for new 
materials in this class is ongoing. With the goal of en- 
larging the space of candidate crystalline structures that 
can potentially support topologically non-trivial insulat- 
ing phases we study in this Letter a class of tight-binding 
models with SO coupling for electrons moving on the py- 
rochlore lattice displayed in Fig. [IJl. Our model belongs 
to the class of 3D 'frustrated hopping' models (TB] and the 
motivation for this study comes in part from our recent 
finding that electrons on the kagome lattice, a canoni- 
cal example of the frustrated structure in 2D, form a 2D 
topological insulator when SO coupling is present fTTJ. 

Our main finding here is that, quite generically, when- 
ever electrons hopping on the pyrochlore lattice acquire 
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FIG. 1: (Color online) a) Pyrochlore lattice is a face-centered 
cubic Bravais lattice with a 4-point basis forming a shaded 
tetrahedron, b) The first Brillouin zone of the FCC lattice 
with high symmetry lines and points indicated, c) Band struc- 
ture of the tight-binding model Eq. |T|. d) Band structure 
with spin-orbit coupling Eq. (|3| for A = —Q.lt (solid line) and 
A = Q.lt (dashed line). 



a band gap from SO interactions the resulting state is 
either a STI or a weak topological insulator (WTI), de- 
fined as a state with j/q = but at least one Z/'i=i,2.3 7^ 0. 
At quarter filling the physics leading to the TI behavior 
on the pyrochlore lattice is somewhat similar to the Fu- 
Kane-Mele (FKM) model on the diamond lattice |T]. SO 
interaction produces Dirac-type spectrum at the three X 
points of the Brillouin zone (BZ), Fig. [T|d, and uniaxial 
crystal distortion is needed to open up a gap. The re- 
sulting Z2 indices are however different from FKM. At 
half filling, the band crossing occurs at the T point and 
is quadratic rather than Dirac-like. In this case SO cou- 
pling by itself opens up a gap and no lattice distortion is 
required. This is, to our knowledge, a unique behavior, 
which produces a highly symmetric 'pristine' STI char- 
acterized by index (1;000). 
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We now supply the technical details supporting these 
claims. Our starting point is the tight-binding model 



Ho 



(1) 



where cj^ creates an electron with spin a on the site 
Tj of the pyrochlore lattice and {ij) denotes nearest 
neighbors. In momentum space Eq. ([ij becomes _ffo = 

with ^fkcr = (cikcr,C2kcr,C3ko-,C4kCT)^ 



and of the form 
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The lower triangle of the matrix is understood to be filled 
so that the matrix is hermitian. The spectrum of 7i^, Fig. 
[ij;, consists of two degenerate flat bands 
two dispersive bands 
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with Ak 



cos (2/03;) cos (2fcy) + cos(2A:a;) cos(2A:2) + 



cos(2fcy) cos(2fcz). E^' touches the two flat bands at the 
r point and the band crossing is quadratic. E^"' and 
E^ ' touch along the lines located at the diagonals of the 
square faces of the BZ. 

At the half-filling, bands 1 and 2 are filled completely, 
and the two degenerate flat bands are empty. This state 
is a gapless band insulator. We now seek terms bi-linear 
in the electron operators that lead to the formation of 
a gap at the quadratic band crossing point. We focus 
on perturbations that do not further break the transla- 
tional symmetry of Hq and preserve T. A natural term 
to consider is a SO interaction of the form 



so 



where A is the SO coupling strength, djj^ are nearest 
neighbor vectors traversed between second neighbors i 
and j, and <t is the vector of Pauli spin matrices. Since 
d^j^ lie in three-dimensional space, the Hamiltonian does 
not decouple for the two spin projections and in k-space 
becomes an 8 x 8 matrix. Fig. [T}i shows the spectrum of 
Hk + Hk*^ . For A > it remains gapless but for A < 
a gap Ago = 24|A| opens at the T point. This peculiar 
behavior can be understood by studying the matrix 7ik + 
H^'^ at k = 0. It is easy to see that the SO couphng splits 
the 6-fold degeneracy into a two-fold degenerate level at 
2t + 16A and a four-fold degenerate level at 2t — 8A. For 
X/t < this allows the four-fold degenerate flat band 
to split off from the two-fold degenerate dispersive band. 
We shall see momentarily that the resulting state is a 
STL 



(3) 



Although the SO interaction reduces the degeneracies 
of bands 1 and 2, they still touch at three inequivalent 
Dirac points = 2'Kf/a, where r = x,y, z. At quarter- 
filling band 1 is fully occupied and it is interesting to ask 
what T-invariant perturbation would open up a gap at 
the Dirac points. We have been able to identify two such 
terms: (i) lattice distortions leading to anisotropy in the 
nearest-neighbor hopping amplitudes, and (ii) modula- 
tions in on-site potentials within the unit cell. Both of 
these preserve the unit cell, the inversion symmetry and 

r. 

For the lattice distortions, since there are six hopping 
amplitudes in the unit cell Fig.[l^, one has many choices. 
We now describe four 'basic' highly symmetric anisotropy 
patterns that open up gaps with equal magnitude at all 
three Dirac points. We then classify the resulting insu- 
lating phases and argue that this classification is in fact 
exhaustive. The basic distortion, labeled by I — 1, 2, 3,4, 
is obtained by selecting site I in the unit cell and changing 
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The -|- sign refers to the six bonds emanating from site 
I whereas the — sign refers to all remaining bonds. This 
can be achieved by deforming the crystal along the axis 
passing through the site / and the center of the tetra- 
hedron. For pattern 1, the Hamiltonian 7i^"^ describing 
this modulation takes the form 
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/ cos{kx — ky) 
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for both spin projections. For / = 2, 3, 4 the signs in front 
of the cosine terms are permuted in an obvious way. The 
full expression for the spectrum of T^k + ^k*"* + ^k"' 
complicated but it is easy to establish that gaps A^is = 
4 1 77 1 simultaneously open up at all the Dirac points. 

As mentioned above a gap also opens up as a result 
of on-site potential modulation. A convenient symmetric 
choice defines pattern / as e; = 3/i and efc^j = — /i with fi 
a constant. 

We now study the topological classes of these insu- 
lating phases. As shown in Ref. [TT] the Z2 topologi- 
cal invariants {vq] ^11^2^3) are easy to evaluate when a 
crystal possesses inversion symmetry. The invariants can 
be determined from knowledge of the parity eigenvalues 
£.2m{^i) of the 2771-th occupied energy band at the 8 T- 
invariant momenta (TRIM) Ti that satisfy = + G. 
The 8 TRIM in our system can be expressed in terms 
of primitive reciprocal lattice vectors as Fi=(„j„2„3) — 
(77ibi + n2b2 -l- n3b3)/2, with Uj = 0,1. Then Va is 
determined by the product (-1)"" = n„j=o,i '^«in2n3' 
and (-l)''i=i.2,3 = 1" -1 '^"-"-"-^ where 6i 
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Dis 


Mass (m^;, mj,, rriz) 




Z2 class 




^2 class 


1 




7? < 


0;100 


r; > 


1;100 


2 


v,-v,v 


ri<0 


0;001 


»7 > 


1;001 


3 


v,v, -V 


r? < 


0;111 


77 > 


1;111 


4 


-V, -V, -V 


r? < 


0;010 


r; > 


1;010 



TABLE L Z2 class for the insulators at quarter filling and the 
corresponding Dirac mass values in the low-energy efi'ective 
Hamiltonian ([5| for different distortion patterns and arbitrary 
A / 0. 



Our model is inversion-symmetric for all perturba- 
tions discussed above and so we can use this method 
to find fs. If we select site 3 of the unit cell, Fig. [l^, 
as the center of inversion then the parity operator acts 
as P['0i(r),V'2(r),V'3(r)>V'4(r)] = [i^ii-r - a-^) , ip^i-r - 
^3)1 V'3(~i'); V'4(~r ~ ^2)] on the four-component elec- 
tron wave function in the unit cell labeled by vector 
r. In momentum space and including spin the par- 
ity operator becomes a diagonal 8x8 matrix Vk = 
diag(e-*'"i '', e'"'^^-^, 1, e-*'^= '')(g)diag(l, 1). It is straight- 
forward to obtain the eigenstates of and the parity 
eigenvalues of the occupied bands numerically, then de- 
termine the Z2 invariants. At half filling, we find that 
(5 = 1 at the T point and S — ~1 at other TRIM, so the 
spin-orbit phase is a (1; 000) strong topological insulator. 

At quarter filling, we find 8 Z2 classes, depending on 
the type of the distortion. 4 of these are STIs and 4 are 
WTIs (Table I). The distinction between STI and WTI 
is easy to understand on physical grounds by considering 
some limiting cases. In the limit rj ^ —t the electrons can 
only move in decoupled parallel planes, each forming a 2D 
kagome lattice. Electrons hopping in the kagome lattice 
with spin-orbit interaction of the form ([3| produce a 2D 
topological insulator [17]. A collection of such planes 
results in WTI even after interplane coupling is restored. 
When 77 — !■ i, on the other hand, the resulting structure 
remains 3-dimensional and STI behavior prevails. 

To develop better understanding for the insulating 
phases at quarter filling we now examine the form of 
the low-energy Hamiltonians governing the excitations 
in the vicinity of the three Dirac points. This is ob- 
tained by linearizing + + H^^^ near and sub- 
sequently projecting onto the subspace associated with 
bands 1 and 2. Near the point we rescale momenta 
as l2Xkx{ky) — > kx{ky) and Atk^ — > fc^, and obtain a 
three dimensional Dirac Hamiltonian, 



2m' 



(5) 



"hC^g are the same with x^y and z permuted in ki and 
<T*. Index I in the mass labels the distortion pattern and 
the values of masses are shown in Table I. We observe 
that I = 1, 2, 3,4 and two possible signs of rj exhaust all 
possible sign combinations for the three Dirac masses. 
Since the Z2 index can only change when at least one of 
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FIG. 2: Band structures for a slab with a III face in various 
insulating phases with A = —Q.lt and \r]\ = 0.2t. The Z2 
index refers to quarter filling. At half filling all four panels 
represent a (I;000) STI. The inset shows the surface BZ with 
high-symmetry points marked. 



the masses goes through zero it follows that our classifi- 
cation in Table I is exhaustive. In particular given any 
cell-periodic pattern of bond distortions and on-site en- 
ergies the Z2 class is uniquely determined by the pattern 
of the Dirac mass signs listed in Table I. 

One can also study the origin of the topologically pro- 
tected surface states using the above low-energy Hamilto- 
nian ([5]). Consider, for the sake of concreteness, a bound- 
ary between two different phases, running along, say, the 
z = plane in real space. We take distortion pattern 2, 
?7 < in the left half-space and pattern 1, 77 > in the 
right half-space. The mass necessarily undergoes a 
sign change across the z — Q boundary. Such a soliton 
mass profile is known to produce massless states in the 
associated Dirac equation, localized near the boundary 
[T^. Specifically, 3D Dirac equation 

[-iT-a, + (a-fc,+a^fc^)T^+m(z)r-]0k(2) = E<t>k{z) (6) 
with m(z —t ±00) = ±mo has gapless solutions 



1 



1 

i 



-/o'™(^')<i2' 



(7) 



extended in the z — Q plane but localized in the trans- 
verse direction, with linearly dispersing energy E± = ±fc, 

[kx 



where (/3k 



iky)/k and k = \ k"^ 



ky. The num- 
ber of gapless states determines the topological class of 
the phases on two sides of the boundary. If the number 
is odd, the boundary is between a WTI and a STI phase. 
If the number is even, the boundary is between two WTI 
or two STI phases. 

To further support our identification of topological 
classes given above, we have performed numerical diago- 
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nalizations of the lattice Hamiltonian Ho+Hso + H^is us- 
ing the slab geometry. Fig. |2] shows the two-dimensional 
band structures of a representative sample of the insulat- 
ing phases obtained for various patterns of crystal distor- 
tions. We consider a slab geometry with two 111 surfaces, 
and plot band energies along lines that connect the four 
surface TRIM. Four bulk bands are clearly visible and 
there are also surface states some of which traverse the 
gap. 

At half fiUing the system is in the (1; 000) STI phase 
(we take A < 0). Irrespective of the lattice distortion 
there is a single Fermi surface (FS) around the T point 
for each surface. Since the two surfaces are inequivalent 
in this geometry the Fermi surfaces are also different. 

At quarter filling more possibilities arise. In the WTI 
phases (0; 111) and (0; 100) the (even) number of surface 
FS depends on the orientation of the surface vector with 
respect to (i^ii^2'^3), as discussed in Ref. |1]. For instance 
the 111 surface has no surface FS in the (0; 111) insula- 
tor while there are two per surface in the (0;100) phase. 
In the STI phase (1; 100) and (1;111) there are 1 and 
3 Dirac points on TRIM, respectively. Near each such 
Dirac point a pair of robust spin-filtered states exists. 
Crossings at other momenta can occur, however there is 
always an even number of such crossings, confirming the 
above arguments. 

Many insulating compounds with pyrochlore structure 
are known to exist [1^]. These follow a formula A2B2OJ 
with A typically a rare earth and B a transition metal 
element. Existing experimental studies so far focused 
mostly on the magnetic pyrochlores due to their promise 
as candidate systems for exotic spin liquid and spin ice 
ground states brought about by the geometric frustration 
inherent to the pyrochlore structure. Our theoretical re- 
sults show that non-magnetic pyrochlores with strong SO 
coupling could exhibit interesting physics. One promis- 
ing candidate is pyrochlore Cd20s207 which shows insu- 
lating behavior below 225K [20j. Band structure calcula- 
tions in this compound favor non-magnetic ground state 
and indicate strong SO effects [21]. Also promising are 
the Ir-based pyrochlores ^2^207 since various Ir-based 
transition metal oxides have been predicted and reported 
to exhibit significant SO effects In addition, 

for A =Nd, Sm, Eu metal-insulator transitions have been 
reported at 36, 117 and 120K, respectively [25] . 

Whether a particular pyrochlore is a topological insu- 
lator can be established only through a detailed band 
structure calculation or an experimental measurement. 
These are clearly beyond the scope of our present study. 
It is however very encouraging to note that, in a preprint 
that appeared after the completion of this work, Pesin 
and Balents [26' derived a semi-realistic tight binding 
model model for A2lr207 {A =Pr,Eu) and found band 
structures for active orbit als closely resembling that dis- 



played in Fig. Id for \/t < 0. They find, in agree- 
ment with our results, that the system at half filling is 
a pristine strong topological insulator. We conclude that 
pyrochlore oxides are likely to open a new frontier in 
the quest for technologically useful topological insulators 
and, more generally, exciting new topological states of 
quantum matter. Clearly, detailed band structure calcu- 
lations and careful experimental studies of these families 
of materials are warranted. 
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